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Abstract : 

We analyze higher twist effect in photoproduction of quarkonium, where the quarko- 
nium is a spin-triplet, S-wave state. We find that the nonperturbative effect of next-to- 
leading twist is contained in three correlation functions related to the initial hadron and 
some of the effect of next-to-leading twist is only suppressed by the inverse of the mass 
square of the quarkonium. A naive estimation indicates that the effect can be significant for 
charmonium. In the analysis we only take the leading order in the small velocity expansion 
for the nonperturbative parts related to the quarkonium. Possible corrections from higher 
orders, especially, corrections from color-octet states are discussed in detail. 

PACS Numbers: 12.38.Bx, 13.85.Ni, 14.40.Gx. 
Keywords: Twist-4 effect, quarkonium production. 



1. Introduction 

Recently the theory of quarkonium production has been developed rapidly(See recent 
reviews given in [1]). Quarkonium production can be regarded as a two-step process, where 
a heavy quark Q and its antiquark Q are produced and then this QQ pair is transmitted 
into a quarkonium. The production of the QQ pair can be treated with perturbative QCD 
because the mass Mq of Q is large, while the transition is a nonperturbative process. 
It is realized that in a rest frame of the quarkonium the quark Q or the antiquark Q 
moves with a small velocity v. Hence a small velocity expansion can be used to treat 
the nonperturbative transition. As a result the effect of the transition can be described 
by quarkonium matrix elements of local operators defined in nonrelativistic QCD [2]. It 
is interesting to note that not only a color-singlet QQ, but also a color-octet QQ can be 
transmitted into a quarkonium. This is in contrast to the color-singlet model and is specially 
of importance to P-wave quarkonium production, because the production rates through a 
color-singlet and through a color-octet are at the same order of v. For 5- wave quarkonium, 
although the contributions from a color-octet QQ is suppressed at least by v'^ in comparison 
with those from a color-singlet, they can be significant [1]. Besides this nonperturbative 
effect related to the transition there is another nonperturbative effect if there are hadrons 
in the initial state. This nonperturbative effect can be analyzed with a twist expansion. 
At the order of leading twist it results in the well-known parton model. The subject of 
the present work is to analyze the effect from the order of next-to-leading twist in inclusive 
photoproduction of quarkonium. 

In the case of inclusive photoproduction of a single jet with a large transversal momen- 
tum one can expect that higher twist effects are suppressed by the power of P^^, since 
Pt is the only large scale in the problem. If one observes instead of a single jet a quarko- 
nium state, there is, in addition to Pt, another large scale-the mass of the heavy quark Mq. 
Naturely one expects that some effects from higher twist may be only suppressed by the 
power of Mq^. In this work we assume that no polarization is observed, so the powers will 
be even numbers. We find that some effect of next-to-leading twist is suppressed indeed by 
Mg^. This indicates that the effect from higher twist may be substantial in the production 
of charmonium, as the mass Mc is not so large. 

In this work we only consider the quarkonium as a S-wave, spin-triplet state and denote 
it as V'. For charmonium and bottonium it corresponds to J/ip and to T respectively. We 
will assume that the nonperturbative parts related to the quarkonium and to the initial 
hadron can be factorized separately and analyze in detail the effect of next-to-leading twist 
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in the i/^-production, where the formation of i/j from a QQ pair is considered at leading order 
of v'^. At this order, only a color-singlet QQ can be transmitted into ip. At higher orders 
of V there are relativistic corrections and contributions from color-octet states of QQ pair. 
We will discuss this in length in a separate section, where we also discuss electroproduction 
briefly. To analyze the nonperturbative effect related to the initial hadron we employ the 
diagram expansion method. This method was first used to classify higher twist effects in 
deeply inelastic scattering[3] , it exactly delivers the same results as those obtained with tra- 
ditional Wilson operator expansion[4] . This method was also applied or extended to other 
processes. An incompleted list of references can be found in [5-10]. However, the extension 
of the method to our case is not quite straightforward. In deeply inelastic scattering, because 
the process is totally inclusive, the perturbative parts of diagrams with different numbers 
of incoming glue as partons can be related to each other with the help of low energy theo- 
rems [3], basically derived from the Ward identity. This leads to that the contributions from 
diagrams with a given number N of incoming glue and from diagrams with N + 1 incoming 
glue can be combined together by replacing a space-time derivative in the corresponding 
nonperturbative part with a gauge covariant derivative. Hence SU{3) gauge symmetry is 
maintained explicitly. The process we consider is seminclusive, where a QQ pair is observed 
indirectly, thus the corresponding low-energy theorems are lacking. Further, the perturba- 
tive parts are much more complicated than those in deeply inelastic scattering. In this work 
we will only consider in the diagram expansion those diagrams whose perturbative parts are 
at the same order of coupling constants and will simply replace space-time derivatives in 
nonperturbative parts by gauge covariant derivatives. 

In our work we will assume that we do not observe any polarization of any particle, and 
we take only leading order in coupling constants. We will also assume that the produced 
ip has moderate transverse momentum, i.e., the t/j is not in the forward region in a fixed 
target experiment. In this region various processes can produce t/j, some of them can not 
be described with perturbative QCD. We neglect heavy fiavor content in the initial hadron, 
which is expected to give small contribution. An interesting aspect of the processes is that at 
this order only glue are responsible for the quarkonium production. Studying these processes 
will help to understand the role played by gluon in hadrons. The process considered here 
was used to study gluon distribution in nucleon before(See [11] and references cited there). 

Our work is organized as the following: In Sect. 2 we introduce notations and the 
diagrams in the diagram expansion which are needed to be considered. The results at the 
leading twist is also given. In Sect. 3 we work out the results for next-to-leading twist. A 
discussion of our results is presented. In Sect. 4 we will discuss possible contributions from 
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higher orders of v, where production of a color-octet QQ pair can also lead to contributions 
for the total production rate. In this section electroproduction is also discussed briefly. We 
summarize our work in Sect. 5. 



2. The Diagram Expansion and Results of Leading Twist 



We consider the inclusive process: 

^{q) + H{P)^^|;{k) + X. (2.1) 

The momenta of particles are given in the brackets. We choose a frame in which the hadron 
H moves in the ^-direction and the photon in the opposite direction. The photon is on-shell 
and H has mass Mh. Throughout our work we will neglect Mh- The correction introduced 
by the target mass may be included by using Nachtmann variable[12]. In oiu' approximation 
the quarkonium mass is twice of Mq. Wc define s = {q + P)"^ and t = (P — k)^. For 
performing the analysis it is convenient to use light-cone coordinate system. In this system 
a 4-vector is expressed as p^' = = {p~^,p~, Pt)- The relation to conventional 

expression is p"*" = -^(p" + P'^)i P~ = '^(P*^ ~ P^)- the system the photon carries the 
momentum = (0,g~,0,0) and H Pi^ = {P~^, ^^,0,0). We introduce in the frame two 
vectors n and £ and a tensor d^^: 

n'^ = (0,l,0,0), £'^ = (1,0,0,0), 

2.2 

df^" = - n'^e" - n^'l". 

Contracting any vector with djf gives its transverse part. We will work in the light-cone 
gauge n ■ G{x) = G'^{x) = 0, where G^{x) = G°''^{x)T°' is the gluon field. In this gauge a 
component of the gluon field strength tensor G^^^ takes a simple form: 

C^+^ix) = -^^G^^^ix). (2.3) 
dx~ 



To analyze the effect upto next-to-leading twist we need to consider the diagrams 
given in Fig.lA, Fig. IB and Fig.lC in the diagram expansion approach. In these dia- 
grams the black boxes represents the nonperturbative part involved by the hadron the 
upper parts, i.e., the blank boxes, contain the perturbative part for production of the QQ 
pair at the leading order of coupling constants and the transition of the pair into t/^, they 
represent the sample of Feynman diagrams given in Fig.2A, Fig.2B and Fig.2C correspond- 
ingly. These diagrams are with Cutkosky cut since we are interested in cross section. With 
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twist-classification, Fig.lA has the leading twist 2, Fig. IB has 3 and Fig.lC has 4. The 
contribution from Fig.lA can be written as: 

74 L r ^4, 



) 



(2.4) 



In Eq.(2.4) corresponds to the upper part in Fig.lA, this part is the contribution from 
diagrams given in Fig.2A. The hadronic matrix element in Eq.(2.4) corresponds to the black 
box in Fig.lA and it is averaged over the spin of H. Because of the color symmetry we can 
write Eq.(2.4) as 

2^'"^ = S'"'^'' -^i^lj fe^-^^^^ 

• J d^xe-'^'-P' < iy(P)|G^'^(a;)G^'^(0)|iy(P) > ^^'^^ 
The contribution from Fig.lC takes the form: 

d^xid^X2d^x^e-'''^-^^-'''^-''^^'''^-^^ 

< H{P)\G'''''''{xi)G''^''''{x2)G'''''''{x3)G''^''''{0)\H{P) > 



(2.6) 



where C is the upper part and corresponds to the diagrams given in Fig.2C. Note C and A 
are at the same order of a and ag- At the leading twist the momenta k'^s in the perturbative 
part C will be set to be proportional to P, so the quarkonium will be produced with zero 
transversal momentum. We will not consider this type of contribution as mentioned in the 
introduction. As for the contribution from Fig. IB, the perturbative part is one order higher 
in Qs than those in Fig.lA. Again we will not consider it. It is worth to point out that 
some contributions from diagrams at higher order of gg will be included in our final results, 
because we later replace space-time derivatives with gauge covariant derivatives for gauge 
invariance. In the case we consider the contribution to the process is only from Fig.lA with 
the perturbative part specified with the Feynman diagrams given in Fig.2A. The main task 
for analyzing twist-4 effect is to separate them from the contribution from Fig.lA. 

To achieve the separation one observes that the space-time dependence of the matrix 
elements in Eq.(2.5) is controlled by different scales in different directions. In — direction 
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this scale is P'^ and it is very large. In +-direction and transverse direction this scale is just 
P~ or Aqcd, the A- parameter of QCD, which are small. With this observation one can 
approximate the x'^- and xx-dependence by Taylor-expanding and xt. This expansion 
is equivalent to the expansion of S^t,{pi) around Pi = Pi = (pi", 0,0,0). The expansion 
reads: 

S,.{p,) = S,Api) + ^^f^\p,^pAP^ 

2c (2-7) 



where Api = pi — pi- The ■ ■ ■ denotes higher order in Api which will lead contributions 
with twist higher than 4. The expansion is called as coUinear expansion. The leading twist 
contribution comes only from the first term in the expansion. This term does not depend 
on p^ and Pit- Taking this term for in Eq.(2.5) we can perform the integration over 
Pi and Pit directly: 

■ ^ / dx-e-'^^^''~ < H{P)\G^^''{x-n)G^-^'{Q)\H{P) > ^^"^^ 

ZTT J 

+ " Higher Twist Terms" 

where we made the substitution p'l = zP'^ . Because of the symmetries of Lorentz boost 
along ^-direction and of rotation around z-axis, the Fourier transformed matrix element can 
be written: 



^ j dx-e-''^^''- < H{P)\G^^^{x-n)G^^^'{{))\H{P) >= 

d^""^ f dx-e-'^^^'^' < H{P)\G^^P{x-n)G\{Q)\H{P) > 

+ n'^n"^ j dx-e-''^^''~ < H{P)\G^^- {x' n)G^^- {Q)\H{P) >, 

(2.9) 

where dj^ is given in Eq.(2.2). In light-cone gauge the dynamically independent gluon fields 
are G^ and G^, G~ can be solved with the equation of motion and the term in Eq.(2.9) with 
G~ will be shown in the next section to be at higher twist. The leading twist contribution 
from Fig.lA can be written now as: 

^'^^^ ^ WT^'^''^^^' ~ I T-^^/^^^^ ■ f"^^^'^ ■ ^^''^^'^^ (2 10) 

+ "Higher Twist Terms" 



5 



The term in [• • •] is just the probabihty for the photon-gluon fusion into '4> and a gluon, 
where the initial gluon is on shell because of = 0. Therefore the interpretation of parton 
model applies. The function fo/niz) can be read from Eq.(2.9): 

fG/H{z) = J dx-e-'^^^-- < H{P)\G'^f{x-n)G\mH{P) > . (2.11) 

Because of the total momentum conservation in + direction the variable z is constrained 
between and 1. The function fa/Hi^) is positive and can be interpreted as the probability 
of finding a gluon in H with the energy fraction z. It is preferred to use strength tensor 
instead gluon field to define fo/Hi^)- Using Eq.(2.3) one can obtain: 

foMz) = J dx-e-'^''^^- < H{P)\G''+^{x-n)G''\mH{P) > (2.12) 

This is the definition of gluon distribution given in [13] in light-cone gauge. For arbitrary 
gauge one needs to supply Wilson line operator between the two tensor operators. 

The perturbative result at the leading twist is well known. For completeness we will 
give this result. For this purpose we introduce 

s={q + Pif, i={k-pif, 

(0) (2) (2.13) 



a^'^^ denotes the contribution at leading twist, u^^-' stands for the contribution at twist- 
4. The invariant s and t by neglecting Mh are simply related to s and t via s = zs and 
I}) = z{t - Ml) 

327r o^o , , ,„ .„.,o C dz 



(i — MV) — z{t — MV) respectively. The result at leading twist is: 



dt 



r d7 



s\s - Mlf + P(t - Mlf + (g + if{s + i - Mlf 



(2.14) 



where Qq is the electric charge of Q in unit e. i?^(0) is the radial wavefunction for at 
the origine, which is related to a matrix element defined in [2]. 

The contribution from the next-to-leading twist is from the term with n^my in Eq.(2.9) 
and the terms with derivatives in Eq.(2.7). We will analyze them in the next section. 
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3. The Results at the Next-to-leading Twist 



We start first with the term with nf^n" in Eq.(2.9). This term corresponds to the 
contribution from the nonpropagating part of the gluon propagator[5] . In this term only 
the — components of gluon fields are involved and they are not dynamically independent. 
With the equation of motion in QCD for these component one can solve them in terms of 
transverse components of gluon fields and quark fields. The equation of motion reads: 

^^^P^ = -^.J«'+(x), (3.1) 

where Dj, = dji' Di, and Dj^ is the covariant derivative in the adjoin representation. 

is the +-component of quark color current. Assuming all fields approach zero at infinite 

space-time, one can solve this equation formally: 

G«'+-(a;+,a;-,XT) = ^- / d^a;- - e){ W^+^l^^, C, XT)r + J"'+(a;+, C, xt)}- 

^ J -oo 

(3.2) 

The function t{z) equals 1 for 2; > and -1 for z <Q. To use this solution we write G~ into 
in the second term in Eq.(2.9) with help of Eq.(2.3) by partial integration. Finally we 
obtain the contribution from this term: 

where the function ei{z) is defined as 

ei{z) = J dx-e--^^-" < HiPMDl^G+^ix-nW + g.J'^'+ix-n)] 

•[(D^G+,(0))« + ^,J«'+(0)]|iy(P)>. 



(3.4) 



This function is positive and has support only for z between and 1. It has a dimension 2 in 
mass while fo/Hi^) is dimensionless. Therefore this term will lead to contribution at order 
of next-to-leading twist. We used subscript Al in Eq.(3.3) to denote this contribution. The 
perturbative part in this contribution is the part in {• • •} in Eq.(3.3) and can be calculated 
straightforward. We obtain: 



d'^Ai 327r |o /nM2 / ei(^) 1 

= ^aa,QQM,\R,m J-'^' ,3(, _ M|)2(t - M|)2(5 + f)^ 

{2s^{s + if - AMl{s^ + sH + 2sP + 2£4) + 2M^(s2 + sP + 4p)} 



(3.5) 



Comparing this contribution with that from the leading twist in Eq.(2.14) there is an 
extra term, the contribution hence may be suppressed by this power. 
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Next we consider the contribution from the terms with derivatives in the colhnear 
expansion in Eq.(2.7). The second term in Eq.(2.7) leads a contribution which can be 



written as: 

j4l r ^4 



• {-%) J d^xe-'^'-P^ < H{P)\{d'^G'''''{x))G''^i'{0) + G''^''{0){d'^G''^''{-x))\H{P) >, 

(3.6) 

where we have used partial integration and translational invariance for the hadronic matrix 
element. Since we do not observe any spin of any particle, by using P(Parity)- and T(Time 
reversal)- symmetries one can show that: 

S^^Api) = S'-^ipi), (3.7) 

while the matrix element is antisymmetric in fiu in contrast to Eq.(3.7). Hence the second 
term in Eq.(2.7) does not contribute. If any polarization is observed, this term will lead to 
a twist-3 effect. 

The third term in Eq.(2.7) leads to a contribution given by: 

• (-1) J dx-e-'^^'Pt < H{P)\G^-''ix)id''d'^G^'{0)f\H{P) > . 



It should be noted that by using the translational invariance one can let the derivatives 
acting on the gluon fields differently than that given above and hence there is an ambiguity 
to form operators in the matrix element. We take the above form since the operator is just 
the trace term subtracted in the corresponding twist-2 operator. The leading contribution 
in Eq.(3.7) comes from those terms in which the contracting Lorentz indices u, a and /3 
are only the transverse indices. If any of them is — , one can show that the contribution is 
at higher twist. For fx or ly being — the argument is the same as for the second term in 
Eq.(2.8). For a or /3 in the matrix element being — one can show because of the covariance 
of Lorentz boost along 2;-direction that the corresponding contribution is suppressed by the 
factor -pzf relative to the terms with a or P being transverse indices. Similarly as to Eq.(2.9) 
one can write the Fourier transformed matrix element in Eq.(3.8) as: 

^ [ dx-e-'^^'Pt <ii-(P)|G'^'^(x)(a"a^G'^'^(0))|iy(P) > 
2n J 

= {d':^d^ + d^^d^ — d'^ d'^^) ■ di{z) ^3 

+ (34"4^ - (4"4^ + <^d^")) • d2{z) 

+ "High Twist Terms" 
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where the functions di{z) and d2{z) are defined as 



^^^^^ = 2^ / dx-e-'^^'^^-lK H{P)\G^^\{x-n){D^^^D'^G\{^)r\H{P) >, 
d2{z) = J dx-e— ^^-"i < H{P)\G'^'+^{x-n){D'j,G\{0)r\H{P) > . 



(3.10) 



In Eq.(3.9) we replaced derivatives with covariant derivatives and dj^dj, with Dj^D^ . The 
symbol {crp} means to take the part symmetric in ap. Calculating moments of the distri- 
butions di{z) and d2{z) one can show that the moments correspond the matrix elements of 
local operators, which are the trace-terms subtracted in the twist-2 gluonic operators. We 
denote the twis-4 contribution from Eq.(3.8) to the cross-section as aA2, which is: 



• {^i(^) • ■ • (^''^ + ^^t"- - d^'df) (3.11) 

ID dp I op I 



Some care must be taken for calculating this part. In S^^{pi) there is a 5-function 
d{{q + pi — /c)^), which is the on-shell condition for the gluon crossing the cut in Fig.2A. 
Therefore the derivatives of S^^{pi) are constrained derivatives. To calculate the derivatives 
we first perform the integration of z with the 5-function in Eq.(3.11), where piT is not zero. 
After the integration is fixed and it depends on piT- Now the derivatives can be easily 
performed, where the indices a and f3 are transversal indices. Finally we rewrite the result 
as an integration of z with the 5-function and set pi = pi. The calculation is tedious but 
straightforward. The results for the perturbative part at order of aa^ are: 



daA2 327r 
dt ~ ~3~ 



aalQQM^\R^ 



dz 



s^{s-MiY{i-MiY{s + i)^ 

d2{z) -, 



(3.12) 
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where the functions /i and /2 are: 



fi{s,i) =8s^£2(s + £)2 + 16Mpt{-s'' - 5sH - AsH'^ + 2sH^ + 6s£^ + U^) 

+ 8M^(s6 + lQsH + 21s'^P - 8sH^ - mH"^ - mi^ - Si^) + 16M^(-s^ 

- 9sH+s^P + 35sH^ + 40s£^ + Ui^) + 8M|(-2s^ + 8sH - 37s^P - 78sP 

- 43f^) + 16M^°(3s3 + UsP + 17F) + 8M^^{-3s^ + 2si - Up), 
/2(s, i) = + 32sH{-2s^ - sH + 2s^P + 2sH^ - 2s£* - 3p) + 16M|s(-2s^ + QsH 

- Ss^t^ - 3st^ - + 32M^(-6s^ - 2sH+s'^t^ + 3s^£^ + 2s£^ + 3£^) 
+ 32M|(4s^ - sH - + sP - 3t^) + 32Ml°(-s3 + sH - sP + P). 



At next-to-leading twist the contribution to the total cross-section for the process in 
Eq.(2.1) is the sum: 



In our results the nonperturbative parts at order of twist-4 in the process are the correlation 
functions 61(2;), di{z) and (^2(2), while that at order of leading twist is fa/Hiz). These 
functions are universal, i.e., they will appear in other processes at certain order of as, 
provided the same hadron H is in the initial state. The functions ei{z), di{z) and d2{z) have 
a dimension 2 in mass and they are at order of ^qcd- '^^^ perturbative parts corresponding 
to these functions are calculated at leading order of a and ctg in this section. Comparing the 
results with that in the last section we can see that the effect from next-to-leading twist is 
suppressed by the inverse of one scale which can be s, t or M^, or any linear combination of 
them. To see more clearly how the twist-4 effect is suppressed we consider two limiting cases, 
these cases can be realized with suitable cuts in experiment. One case is the kinematical 
region where s, |^| >> M^, but t is not so large that the 4) does not appear in the forward 
region, whose definition will be given later. Another is for s >> 

Case 1). s, \t\ » M|. In this case, one can neglect any power of and In this 
approximation the expression for the differential cross-section simplifies and the contribution 
related to d2{z) can be neglected: 



(3.13) 



(3.14) 




(3.15) 



+ + {s + if 

s2(s + £)2 



•/g/h(^) + 2 




+ 8 



di{z) 



}■ 
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From Eq.(3.15) we clearly see that the twist-4 effect is suppressed only by , but not by 
the inverse of any large scale of s or i. 



Case 2). s >> M|, In this limit we neglect any power of and |. The differential 
cross-section becomes 



da 327r o ^ ^ , ^ ,^^,9 f dz 1 
-aalQlM^\R^{Q)\^ j — 



dt 3 ^' 7 z s^{i-M^)^ 



(3.16) 



Again the twist-4 effect here is suppressed by ^ or by the inverse of a linear conbination 
of M| and i. 

In general there is always some higher twist effect which is suppressed only by . 
From the discussions of the above two cases one can not make this effect negligible through 
suitable cuts in experiment. Therefore the twist-4 effect may be substantial for charmonium 
production, as the mass Mg is not so large. Because the nonperturbative functions 61(2;), 
di{z) and d2{z) are unknown, one does not know how large the twist-4 effect can be. A naive 
estimation could be given by noting the fact that the effect we considered is basically from 
the correction of the parton-model, where one takes the off-shellness and the transversal 
momentum of the parton into account. With this fact one may make an Ansatz for a naive 
estimation as follows: 

61(2;) = d^iz) = d2{z) = AVg/h(-2), (3.17) 

where the parameter A can be regarded as an average of the transversal momentum of the 
initial gluon and it should be at the order corresponding the radius of the initial hadron. 
With this Ansatz the differential cross-section upto twist-4 can be written: 



da 327r n 
-dt = ^""^ 



/dz 
— {W2{s,t) + A^W4{s,t)) ■ fG/H{z). (3.18) 



We will give some numerical results for the quantity p defined as 

_ W2{s,i) + A'^W4{s,i) 
W2{s,i) 



(3.19) 



If the twist-4 effect is negligible, p should be close to 1. An useful variable used commonly 
in experiment is 2, it is defined as: 

h ■ P 

z = (3.20) 
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which should not be confused with the intergration variable z in previous equations. If 
is produced with z > 0.9, it is in the froward region. We plot the quantity p as a function 
of z for a given s, where we take A = SOOMeV, this value may correspond to a proton as 
the initial hadron. The results are given in Fig.3 for s = (lOGeV)^ and for s = (30GeV)^ 
respectively and the quarkonium is J/V'- On can take these as representative estimations 
for ^/s = 30GeV and ^/s = lOOGeV. From the figure we note that the twist-4 effect becomes 
smaller when s becomes larger and it is more significant in lower 2;-region than that in higher 
2;-region. For s = (SOGeV)^ it can be an effect at 30%. 

4. Higher Order Effect in the Small Velocity Expansion 

In previous sections we analyzed the nonperturbative effect at twist-4 related to the 
initial hadron, while the nonperturbative effect related to the quarkonium is only considered 
at leading order of the small velocity expansion. In this section we discuss higher order effect 
in this expansion. At leading order of v only a color-singlet QQ pair in a ^Si state can be 
transmitted into the ^Si quarkonium, the probability for this can be regarded as |i?^(0)p, 
which should be related to a matrix element defined in NRQCD[2]. The relation is 

l^^(0)P = y <0|Ofe5i)|0>. (4.1) 

We will not give the definition of the matrix element and of those appearing later. For the 
definitions we refer reader to [2]. We will take the leading order here as Corrections to 
the leading order results will start at order of v"^. 

At next-to-leading order, i.e., at v'^, there is only the relativistic correction, which takes 
the effect of the relative motion of Q and Q into account. This effect is analyzed within the 
color-singlet model in [14], where the correction is parameterized in term of the parameter 
e defined as: 

= 2Mq + e. (4.2) 

Numerical results given in [14] show that for charmonium the correction at = 14.7GeV 
is significant in the high z region, where z is defined in Eq.(3.20) and e = 0.16Mc. However, 
there is a uncertainty in choosing the parameter e. If one thinks that the binding force is 
Couloumb force, e defined in Eq.(4.2) is negative. One can also show by starting the energy- 
momentum tensor the relation that e equals to the negative matrix element in NRQCD of 
the kinetic-energy operator of Q or Q inside the quarkonium, if Mq in Eq.(4.2) is the pole- 
mass[15]. However, this relation holds at the tree-level and at the leading order of v. All 
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these show that e can also be negative. With a negative e the correction turns into negative 
too. It should be pointed out that the analyze in [14] is done with color-singlet model. 
With the new NQRCD factorization in [2] the nonrelativistic correction is parameterized 
with other type of matrix elements instead of e and a re-analyze seems needed. 

At order of v^, except the relativistic corrections there are also other corrections due 
to the fact that a quarkonium state consists many components, where a QQ pair carrying 
the quantum numbers of the quarkonium is the dominant one. Other components consists 
fo a QQ pair and soft glue, the probability to find them in a quarkonium is suppressed by 
powers of v relatively to that of the dominant one. Because of this a produced QQ with 
quantum numbers other than those of the quarkonium can become one of these components 
by combining soft glue and hence the production of such QQ pair will also contribute to the 
production rate of the quarkonium. In our case such QQ pair can be a color-octet QQ pair 
in ^Sq, ^Si or ^Pj state, where J = 0, 1,2. We use the notation QQ(8)(^^~^^Lj) to denote 
these states. With the velocity power counting rule given in [16] and the method in [2] one 
can show that the contributions from these color-octet QQ pairs are at order of u^. Taking 
the state QQ(8){^Pi) as an example, the probability to find such state combined with a soft 
gluon as a component of the quarkonium is at order of v^, while the formation of Q and Q 
into a P-wave state has a chance at order of v^, hence the production of the quarkonium 
through such state will happen at order of v'^. 

We first discuss the contribution from a QQ{8){^Si) state. Such state can be produced 
not only by photon-gluon fusion at order aa^, but also by photon-quark scattering at the 
same order. However, numerical results show that the contribution from the photon-quark 
scattering is negligible. It is interesting to note that the contribution from the photon-gluon 
fusion is represented by the same diagrams for the color-singlet contribution worked out 
before. To include this contribution one needs only to replace the matrix element in our 
previous results via: 

< 0\OtCSi)\0 >^< O|0f (35i)|0 > < 0\Ot{^Si)\0 > . (4.3) 

where the NRQCD matrix element < 0|Og ("^S*!) |0 > is for the transition of the color-octet 
QQ. Due to the uncertainties in matrix elements and others like effects in higher order of 
as, i^-factor problem • • • etc, it is hard to identify this color-octet contribution. 

The contributions from the QQ(^8){^So), QQ(8){^Pj) state at leading twist are analyzed 
in [17,18,19] for J/ip. The matrix element for the transition of these states into J/^, 
i.e., < 0\Qt{^So)\0 > and < 0\Of{^Pj)\0 >, is estimated[19,20]. The situation here is 
controversial. In the forward region, the contributions are from the subprocess — > Q+Q 
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at the order of (xcXg. In experiment the measured cross section contains also contributions 
from elastic diffr active processes which can not be described with perturbative QCD. That 
means that the cross section estimated from the subprocess should be smaller than the 
measured one, if one neglects all higher order effects. A comparison with experiment for 
J/ip production shows that the predicted cross section from the subprocess [17] and also 
including color-singlet contributions[18] is much larger than measured one, if one uses the 
numerical values for matrix elements estimated from experiment at Tevatron[20], they are 
roughly: 

<0\QtCSo)\0> 

Ml = 2JTT '° 

It seems that these matrix elements arc overestimated at least by one order of magnitude. 
Various reasons for that discrepancy are proposed like corrections from higher order of as 
and uncertainties in extracting numerical values of the matrix elements. One should also 
note that the predictions may suffer large uncertainties because of twist-4 effect, as this effect 
is only suppressed by in the forward region. In addition, the relativistic correction at 

order of v"^ discussed before is significant in the froward region and can be negative. 

In the nonforward region, these color-octet states can be produced at order of aa^. 
At leading twist the production is through photon-gluon fusion described by Fig.lA, where 
the blank part represents not only the Feynman diagrams given in Fig.2A, but also those 
like the one given in Fig.4, in which the gluon in the final state is emitted by the incoming 
gluon. A comparison with available experiments is made in this region in [17,18]. It is 
shown that the contribution from these color-octet states becomes dominant for 2; > 0.5 and 
including this contribution leads to that the differential cross section increases rapidly as 
z increases. This rapid increasing however is not favoured by experimental data[17,18]. If 
the relativistic correction discussed before is negative, the increasing may become mild. An 
important question related to our work is whether there is some twist-4 effect suppressed 
only by Mq as in the color-singlet contribution or not. The answer is no. To see this 
one needs to understand why the twist-4 effect analyzed before is suppressed by Mg^. If 
one expands the perturbative part at leading twist in Mq, the leading power is 2. In the 
final result in Eq.(2.14) an extra Mq^ comes in because of the normalization. The leading 
power of the perturbative parts at twist-4 is 0. Therefore there is twist-4 effect suppressed 
only by Mg^ in the production of a quarkonium. For the above color-octet states the 
leading power of the perturbative parts at leading twist is zero or negative, and in calculating 
perturbative parts at twist-4 extra negative power in Mq can not be generated, so there 
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will be no twist-4 effect which is suppressed only by Mq^ like those in Eq.(3.15). 

Before ending this section we briefly discuss electroproduction. This is basically the 
process 'y* + H ^ tj; + X where the photon 7* is off-shell and the off-shellness is —Q^. One 
can perform an analysis of twist-4 effect as we did here for a real photon. One can show 
for a fixed Q"^ that there is no twits-4 effect suppressed by Mq^. The reason is the same 
as for the color-octet states mentioned above. In photoproduction the forward region is the 
kinematical region where the produced quarkonium moves roughly in the same direction 
as the photon beam in a fixed target experiment. In this region the contribution from 
elastic and diflFractive processes to the production rate is significant. In electroproduction 
like HERA experiment the produced quarkonium in the forward region can have large angle 
in respect to the hadron beam, one may expect that the contribution from elastic and 
diflFractive processes may not be significant. If the events in this region can be described by 
perturbative QCD, the theoretical predictions indicate that there is a peak in the forward 
region, resulted by color-octet states produced by the subprocess 'y* + G ^ Q + Q and 
the subprocess at twist-4 described by Fig.lC. In experiment NMC data of J/V'-production 
shows that there really is a peak in the forward region[21]. However, the data was given as 
the rapidity distribution, hence many information was lost for a detail study of the peak. 
On the other hand, it is also pointed out in the second reference in [11] that the peak can 
be well explained by the theoretical results at leading twist. If more data in experiment is 
available, this region may be a good place to study the effect from those color-octet states 
and from twist-4. 

5. Summary 

In this work we have analyzed twist-4 effect in photoproduction of a ^Si quarkonium. 
The nonperturbative parts at order of twist-4 are factorized into 3 correlation functions 
related to the initial hadron. The perturbative parts are calculated at leading order of 
coupling constants. We find that the twist-4 eflFect suppressed by Mq^ in comparison with 
that at leading twist. Hence the twist-4 effect in photoproduction of J/ip may be substantial, 
and it can not be made negligible by experimental cuts. Naive estimations indicate that the 
effect is significant at ^/s around several ten GeVs. However, a definite conclusion requires 
a detailed study. 

In the analysis we only take the leading order of the small velocity for the nonpertur- 
bative eflFect related to the quarkonium. The eflfects from higher orders of v upto are 
discussed in detail. Especially, at order of production of various color-octet states of QQ 
can contribute to the total production rate. However, the contributions from these states 
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will not receive from twist-4 corrections suppressed only by Mq^ in the nonforward region, 
except the contribution from the color-octet state QQ{8){^Si). But it is hard to distinguish 
this contribution from those of the color-singlet as discussed in the last section. We have 
also discussed some aspects in electroproduction, here the twist-4 corrections are also not 
suppressed by Mq^. It seems that the appearance of such suppression in photoproduction 
of ^Si quarkonium may be a special case. 

Unlike the gluon distribution fa/Hi^)-! which is measured in experiment in many places, 
the three correlation functions for the nonperturbative effect of twist-4 are unknown and 
need to be measured. These three functions are defined in Sect. 3 in the light-cone gauge. In 
arbitrary gauge one needs to supply Wilson-line operators between the gluonic operators in 
these definition. It should be noted that the definitions in Sect. 3 are unrenormalized versions, 
i.e., at tree level. A definition in renormalized version like that for correlation functions at 
twist-2 given in [13] may be difficult. The main obstacles are that these functions may have 
power divergences and the perturbative parts corresponding to correlation functions at twist- 
2 may have problem with convergence due to possible presence of renormalon. The later 
indicates that there may be still some nonperturbative effect contained in the perturbative 
part at leading twist. A perfect definition would be given so that the two problems would 
be solved(See discussions in [22]). 

Finally, it should be noted that the twist-4 effect in the kinematic region considered 
here is basically from the corrections of the parton model, where one takes the transverse 
momentum and the off-shellness of incoming partons into account. With this in mind it 
is possible to make a model for a more realistic estimation of twist-4 effect than the naive 
estimations given in Sect. 3. Currently a model is under development in hope to study the 
effect in detail. 
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Figure Caption 



Fig.lA: The diagram appearing in the diagram expansion. This diagram leads to con- 
tributions at leading twist and at higher twist. The wavy lines are for gluons. The lines 
connecting the black box with the upper part denote only momentum flows and contractions 
of color- and Lorentz-indices. 

Fig. IB: One of the two diagrams leads to contributions at next-to-leading twist. In these 
diagrams the upper part is connected to the black box with three gluon lines. 

Fig.lC: The diagram leads to contributions at next-to- leading twist, the momentum flows 
are given with k4 = ki + k2 — k^. 

Fig.2A: A diagram corresponds to the upper part in Fig.lA. Other diagrams are obtained 
through permutations of gluon- and photon-lines. The line is for heavy quark. The black dot 
denotes the nonperturbative transition into a quarkonium. The thick line is for quarkonium. 

Fig.2B: A diagram corresponds to the upper part in Fig. IB. Other diagrams are obtained 
through permutations of gluon- and photon-lines. 

Fig.2C: A diagram corresponds to the upper part in Fig.lC. Other diagrams are obtained 
through permutations of gluon- and photon-lines. 

Fig.3: The quantity p as function of z deflned in Eq.(3.20). The line is for s — (lO)^GeV, 
the dashed one is for § = (30)^GeV. 

Fig. 4: One of possible diagrams in addition to those in Fig.2A contributes to the production 
rate through color-octet states 
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